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ABSTRACT 

Recent interferometric observations of the young stellar object(YSO) HH 30 have 
revealed a low velocity outflow in the 12 CO J=l-2 molecule line (Pety et al. 2006). We 
present here two models of the low velocity disc winds with the aim of investigating 
an origin of this molecular outflow. Following Andlada et al. (2006) we treated HH 30 
as a binary system. Two cases have been considered: i) the orbital period P = 53 yrs 
and ii) P ^ 1 yr. Calculations showed that in the first case the outflow cone had a 
spiral-like structure due to summing the velocities of the orbital motion and the disc 
wind. Such a structure contradicts the observations. In the second case, the outflow 
cone demonstrates a symmetry relatively to the system axis and agrees well with the 
observations. 
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1 INTRODUCTION 

The Herbig-Haro object HH 30 is a well-known YSO with a 
flared circumstellar (CS) disc seen practically edge-on (Bur- 
rows et al. 1996; Ray et al. 1996) and highly collimated op- 
tical jets propagating perpendicular to the disc, firstly ob- 
served by Mundt & Frid (1983). The star itself is not visible, 
it is detected only through radiation scattered by the CS 
disc. According to Pety et al. (2006) (hereafter P06), HH 
30 is a star with a mass of 0A5Mq, spectral type around 
of Ml and age of (1 - 4)10 6 years. The millimeter observa- 
tions showed that it is surrounded with the CS Keplerian 
disc with the radius of about 420 AU and the mass of about 
4 ■ 1O~ 3 A40. Optical observations showed the disk of about 
450 AU in diameter (Burrows et al. 1996). 

Since 1996 this object has been a subject of intense 
investigations by different authors (Bacciotti et al. 1999; 
Stapelfeldt et al. 1999; Wood et al. 2000, 2002; Cotera et al. 
2001; Watson & Stapelfeldt 2004). Optical observations with 
the Hubble Space Telescope (HST) as well as the millimeter 
interferometry revealed a strong variability and asymmetry 
in the radiation coming from the CS disc. Burrows et al. 
(1996) detected changes in the brightness ratio of the upper 
and lower nebulae between 1994 and 1995, and the presence 
of a weak lateral brightness asymmetry. Stapelfeldt et al. 
(1999) found that in 1998 the right (north-northwest) side 
of the upper reflection nebula brightened dramatically, and 
the other side faded relative to 1995. The system had ap- 
peared nearly symmetrical about the jet axis in 1995, but 



E-mail: tamb@gao.spb.ru (LVT); grinin@gao.spb.ru (VPG) 



in 1998 one side of the disc was four times brighter than the 
other. 

The interpretation of the HH 30 brightness variability 
and asymmetry has a great importance. Stapelfeldt et al. 
(1999) excluded large-scale motions of the matter in the 
outer disc as a possible reason of such variations. They ar- 
gued that the outer disc acts as a screen, on which moving 
illumination patterns are projected from the inner disc or 
the central star. They suggested two mechanisms which are 
able to form such light changes: (i) illumination by the star 
with bright accretion hot spots (Wood & Whitney (1998)) 
and (ii) presence of voids or clumps in the inner disc casting 
beams or shadows onto the outer disc. Each of these mech- 
anisms can produce asymmetry in the nebulae brightness, 
which shifts from one to another side of the disc. Watson & 
Stapelfeldt (2004, 2007) did not find the period to be cor- 
related with the stellar rotation period or any other period. 
According to them, the variability of HH 30 has the more 
complex character and caused by the changes in the CS ex- 
tinction in the inner disc. Tambovtseva and Grinin (2008) 
showed that such extinction variations can be caused by the 
dust component of the azimuthally non- homogeneous disc 
wind. A knowledge of the asymmetry mechanism would give 
a possibility to determine in detail a structure of the nearest 
regions of the young star that has not been done with the 
help of the telescopes. Also this permits us to study more 
definitely the physical processes leading to formation of such 
structures. 

One of the main result of the HH 30 study obtained 
during recent years, was the detection of the slow matter 
outflow in the 12 CO(2-l) molecule line by P06. This outflow 
was highly asymmetric (one-sided) in spite of the presence 
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of relatively symmetric jets, and originated from the inner 
parts of the disk, had a conical morphology with a 30° half 
opening angle and a radial velocity of about 12 km s _1 . 
P06 concluded that the observed CO emission arose from 
material that has been directly launched from the disc and 
moved ballistically in the first few 100 AU from the star; the 
re-collimation could occur at the distances larger than 1000 
AU. According to their interpretation, the high-velocity op- 
tical jet and low-velocity outflow are associated with differ- 
ent parts of the disc wind. The former corresponds to its 
densest part (< 1AU) and the latter to its outer part (5 - 15 
AU). The mass loss rate of the outflow M m was estimated as 
6.3 x 1O _8 M0 yr _1 ; it is almost 60 times more than that of 
the jet equaled to fa 1O _9 M0 yr~ x . Because however of their 
velocity differences it follows that a majority of the kinetic 
energy of the disc wind (L w fa 6 ■ 10 31 erg c _1 ) is contained 
in the narrow collimated jet, while the bulk of the mass loss 
occurs via the low-velocity wind component. 

Recently Anglada et al. (2007) revealed periodic 
changes in the jet's trajectory and connected this with the 
binarity of HH 30. They concluded that the structure of 
the HH 30 jet can be described by a wiggling ballistic jet, 
arising either by the orbital motion of the jet source (i.e., 
the secondary) around the primary or by precession of the 
jet axis because of the tidal effects of the secondary. In the 
first scenario, the orbital period is 53 years, and the sum of 
masses is 0.25 - 2 M . In the second scenario, the mass of 
the jet source (i.e., the primary) has to be ~ 0.1 — IMq, 
the orbital period is less than 1 year, and the mass of the 
companion less than a few xO.O1M . In both scenarios the 
separation of the binary does not exceed 18 AU. Knowing 
the size of the flared disc the authors conclude that this is 
circumbinary rather than circumstellar disc. 

In this paper we model the disc wind in the binary 
system in the framework of the two scenarios suggested by 
Anglada et al. (2007). In section 2 we briefly formulate the 
problem for the first scenario and present results of calcula- 
tions, in section 3 we do the same for the precession scenario. 
Then we discus results and make the conclusion. 



2 DISC WIND IN THE YOUNG BINARY 
SYSTEM 

We consider the wind from the low-mass secondary com- 
panion moving in a circular orbit in the young binary sys- 
tem. The companion accretes matter from the remnants of 
the protostellar cloud, the so-called circumbinary (CB) disc. 
In detail the model is presented in our early paper (Grinin 
& Tambovtseva (2002)). Here is a brief description of the 
model. Following Artymowicz & Lubow (1996), we treat 
the low-mass companion as the main accretor. The wind 
particles, consisting of gas and dust in the standard propor- 
tion, are ejected from the surface of its accretion disc within 
a cone having fixed inner and outer angles. It is supposed 
that in the frame of the secondary, the outflow has an az- 
imuthal symmetry. In the coordinate system of the primary, 
the wind becomes asymmetric due to the vector summing 
the particles and the orbital motion velocities. Trajectories 
of the wind particles in the gravitational field of the primary 
have been calculated in the ballistic approach. We imitated 
the quasi-continuous process by ejecting the particles in a 
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Figure 1. The distribution of the wind particles in the com- 
mon (with the primary) envelope created by the disc wind of the 
secondary component in the models A, B and C (the XZ pro- 
jection). The left panel: the computed particles distribution; the 
right panel: the same after smoothing procedure (details are in 
the text). Letters refer to the model. Distances X and Z are given 
in units of the orbital radius. 

small and equal time interval. Computation stopped when 
the particles reached the outer radius R ou t or the equatorial 
plane of the binary. 

In order to reproduce the molecular outflow in HH 30 
we used parameters estimated by P06: the outflow matter 
is mainly located on the thin walls of a cone with an semi- 
opening angle of 30°. The radial velocity V w = 12± 2 km/s, 
the tangential velocity component is absent. The mass loss 
rate M w = 6.3 x 1O~ 8 M yr" 1 . 

Taking this into account, we adopt the following param- 
eters for our model: mass of the primary M* = O.45M , and 
the eccentricity e = 0. In order to provide the orbital period 
of 53 yrs, the radius of the orbit vq = 11 AU has been cho- 
sen. The initial radial velocity of the wind particle, in units 
of the corresponding Keplerian velocity, was an input pa- 
rameter and varied from 2 to 5. The Keplerian velocity at 
r = 11 AU is about 6 km s _1 . Three models of the outflow 
are presented in Fig. 1, and their parameters given in the 
Table (models A, B, C). 

Thus, in the models considered below, the common pa- 
rameters are the parameters of the orbit and a zero tangen- 
tial component of the wind velocity U. The models differ 
with radial velocities of the wind particles V and angles of 
their ejection 9 (Table 1). The angles are counted from the 
symmetry axis. Figure 1 (the left panel) displays the com- 
puted distribution of the wind particles in the common en- 
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Table 1. Parameters of the models 



Model 


e 


V 


U 


A 


30 


3 





B 


45 


3 





C 


30 


5 





D 


30 - 40 


2 


1 
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Figure 2. Isodcnsity contours for models A, B and C. The details 
arc in the text. 



velope in projection on the XZ plane. The right panel of Fig. 
1 gives the same distribution after the smoothing procedure, 
which has been made as follows. We set a "theoretical beam" 
with a spatial resolution equaled to that used by P06, when 
observing the emission in the 12 CO(2-l) molecule. Then, us- 
ing this beam we scanned a theoretical image shown on the 
left side of Fig. 1. We show only upper planes in the figure. 

An increase in the particle ejection angle (models A and 
B) leads to more amorphous structures, while an increase in 
the initial radial velocity of the wind (models A and C) to 
the larger step of the spiral, but a common structure per- 
sists. The radial velocities estimated from the observations 
do not permit us to increase the start velocities infinitely. 
We checked also the other variants (low radial velocities, 
random ejection angles), but a cone structure has not been 
obtained in the framework of this scenario. 

If we assume, as has been done by P06, that the emis- 
sivity of the gas is proportional to the local density, and 
the brightness distribution is proportional to the column 
density, one may consider the isodensity map as an inten- 
sity map. Such maps are given in Fig. 2 for models A, B 
and C. Contours of equal densities were calculated using 
the smoothed distributions of the prob particles. The con- 
tours are shown in the grey scale where the outermost curve 
corresponds to the level of 100 probe particles, and the max- 
imal level is about 3000 - 4000 particles depending on the 
model. Each level differs from the next one of 100 probe par- 
ticles. The method of the transition from the probe particles 
to the real column densities (in cm~ 2 ) is described in our 
paper (Grinin & Tambovtseva 2002). 

Comparison of the isophotes obtained for these models 
with the observed ones (Fig. 1 in the paper by P06) leads to 
the conclusion that this scenario is not realized in the HH 
30 system. 




Figure 3. The wind geometry. 



3 PRECESSION SCENARIO 

Anglada et al. (2007) suggested a second scenario: the wig- 
gling jet structure is caused by precession of the jet axis due 
to tidal perturbations from the low mass companion. In this 
case the binary system is very close. The jet source is an 
accretion disc of the primary and the question arises: where 
does the low velocity wind come from? P06 concluded that 
for a single star, the source of the biconical molecular outflow 
lies on the surface of the accretion disc at 5 - 15 AU from the 
center. This means, that in the close binary with the low- 
mass companion, located at the distance of < 1 AU from the 
central star, the inner disc of the primary cannot be a source 
of the molecular outflow, because the launched region would 
be near the inner boundary of the CB-disc. In the case of 
a circular orbit, one can assume that the inner boundary 
of the CB-disc possess an axial symmetry. Therefore, the 
initial conditions for the disc wind will not depend on the 
azimuth, and, thus, an outflow produced by such a wind 
will be axially symmetric. In this scenario, the wind parti- 
cles simultaneously start from the inner CB-disc in small 
time intervals. As above, calculations have been made using 
a ballistic approach. Trajectories of the particles have been 
calculated exactly as in the first case, differing only in initial 
and boundary conditions, and with that the orbital motion 
of companions did not influence anyhow the particles kine- 
matics. The radius of the orbit ro = 0.75 AU. With such a 
radius and the mass of the star 0.45 Mq, the period of the 
system is about 1 year. According to calculations by Arty- 
mowicz & Lubow (1994) the inner radius of the CB-disc res 
is 3-4 times greater than the orbital radius ro. We adopted 
res = 4ro ~ 3AU. In this case the Keplerian velocity at the 
inner radius of the CB - disc is about 13 km/s. 

Table 1 includes parameters of the model D considered 
in the second scenario. The radial and tangential compo- 
nents of the initial velocity of the wind particles are ex- 
pressed in units of the Keplerian velocity at res- The par- 
ticles are ejected in the cone constrained with the angles 
6 1 — 40° and 62 = 30° (Fig 3). As calculation show, in this 
case the disc wind produces axially symmetric outflow with 
parameters close to those observed. The isodensity map cal- 
culated for such an outflow is presented in Fig. 4. 

Let us estimate the optical depth of the disc wind r 
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Figure 4. The same as in Fig. 2 for the model D. 



in this model of HH 30. Denote A9 = 9i — 9i- Taking into 
account, that the thickness of the molecular outflow cone is 
small (P06), we assume that 9i « 9i = 9. Then, A6> < 9. 
With the help of the Eq. (1) from the paper by Grinin & 
Tambovtseva (2006), one can obtain the optical depth of the 
cone "wall" along the radius which at the small 7 can 
be replaced by the radius r\Ti which is parallel to the disc 
plane and intersects the disc axis at the distance h from the 
point O' (Fig. 3): 

_ M w K\ . . 

T ~~ "K7 2-Khsm9 ' 1 ' 

Here M w and V w are the mass loss rate and the velocity 
of the disc wind material respectively, n\ is a coefficient of 
absorbtion by the dust component of the disc wind at the 
wavelength A. It should be noted, that when deriving this 
formula we took into consideration, that the outflow in HH 
30 occurs with a strong deviation from mirror symmetry: 
the major part of the matter flows in the upper cone. 

Substituting the values of M w and V w in the expres- 
sion (1) and adopting 9 = 30° and k\ — 250 cm 2 /g (that 
corresponds to the wavelength ~ 5000 A for the CS dust), 
we obtain the optical depth in the base of the wind (h = 5 
AU) r sa 3. Thus, in the model of the outflow adopted in 
this section, the optical depth for the base turns out to be 
of the order of few units. We estimated the angle 7 with the 
help of the scale height law for HH 30 (see Table 3 given 
by P06). We obtain that the peripheral parts of the flared 
disc are illuminated by the stellar radiation under the angle 
7 = 10 — 15 degrees. This means, that at the higher regions 
(the radius rjrj in Fig. 3), the disc wind could still be opaque 
for radiation from the star. If the wind has an azimuthally 
non-homogeneous structure, this may result in the irregular 
illumination of the disc, and thus in the variability of the 
radiation scattered by the disc. 

We also checked the effects of such parameters as radial 
and tangential velocity components on the wind density dis- 
tribution, and did not present these results since all of them 
were analogous to that discussed above. An increase in the 
initial radial velocity gives a more extended outflow, the 
greater tangential component of the initial velocity leads to 
the greater opening angle, but a principal structure remains 
the same: there is a concentration of the wind particles to- 
wards the cone's walls. 



4 DISCUSSION 

Let us put aside the problem of the binarity of HH 30 and 
discuss other possible mechanisms for the formation of the 
molecular outflow. A comprehensive overview of the molec- 
ular outflow models is presented by Arce et al. (2007) (see 
also the references therein). All of them can be separated 
into four broad classes: (1) wind-driven shells, (2) jet-driven 
bow shocks, (3) jet-driven turbulent flows, and (4) circula- 
tion flows. In the first three models the molecular outflow is 
entrained by an underlying wind or jet, in the fourth, it is 
formed by deflected infalling matter. 

Pety et al. (2006) discuss the second entrainment sce- 
nario as a possible mechanism of the outflow origin. Ac- 
cording to this mechanism, the observed outflows consist in 
ambient molecular gas that has been put into motion by 
large bow shocks propagating down the underlying jets (e.g. 
Gueth & Guilloteau 1999). This mechanism explains well 
outflows observed on large scales (3> 1000 AU). In the case 
of HH 30, the outflowing molecular gas is continuously col- 
limated down to the very close vicinity of the star, on spa- 
tial scales that are smaller than the disk size. This implies 
that the outflow arises from material that has been directly 
launched from the disk surface. Also, it is difficult to under- 
stand why CO emission is absent in the southern lobe. In 
the framework of the jet-driven flow models this means that 
properties of the ambient medium in the north and south di- 
rections differ strongly . Thus, the main arguments against 
entrainment and pro disc wind scenario are a conical geome- 
try, a radial velocity near the launching region, and a lack of 
CO emission in the south lobe. For these reasons, P06 ruled 
out the jet-driven flows from considerations in favor of disc 
wind mechanism. Q. 

Nevertheless, if we assume that the molecular outflow 
originates from the outer part of the disk wind, launched at 
a few AU from the star, another problem arises: how to ex- 
plain the interval between the jet and the outflow, in other 
words, the lack of detected material at intermediate veloci- 
ties and angles? The binary status of HH 30 helps to avoid 
this difficulty. In both models of the binary system consid- 
ered in the present paper, the jet and the disc wind originate 
from different sources. It should also be taken into account 
that in the binary system with the low-mass companion a 
cavity free of matter arises due to the gravitational pertur- 
bations (Artymowicz & Lubow 1994), therefore, this pause 
may be naturally explained. 

Based on the analysis of the observed (Hartigan et al. 
1995) and computed forbidden lines, Kwan & Tademaru 
(1995) conclude that a substantial fraction of the emitting 
wind particles has low velocities. Can the low-velocity part 
of the disc wind, originating from the warm disc (Ferreira, 
Dougados & Cabrit 2006, Ferreira 2007), produce such a 
molecular outflow as observed in HH 30? This question needs 
to be explored. In the case of a single star a discontinuity 
in launch regions of jets and disc winds would be a more 
serious issue than in the case of the binary system, but the 
disc wind problem in binaries remains to be investigated. 



1 It should be noted that entrainment mechanism would lead 
to the conical structure of the molecular outflow but arguments 
mentioned above prevent its application. 



Disc wind in the HH 30 binary models 5 



5 CONCLUSION 

We calculated the trajectories of the disc wind motion using 
a ballistic approach within the framework of the two binary 
system scenarios for HH 30. In the first model, the disc wind 
originates from the accretion disc of the secondary compan- 
ion; in the second model, the source of the low-velocity disc 
wind is the inner region of the circumbinary disc. 

Calculations showed that in the first case, the matter 
outflow had a spiral-like structure due to the orbital motion 
of the wind source around the primary. Such a structure 
does not agree with the axially symmetric cone morphology 
of the outflow in HH 30 observed by Pety et al. (2006). The 
cone outflow geometry is typical for the second model, in 
which the well collimated and highly accelerated jet, and 
the slow molecular outflow, originate from the different re- 
gions of the binary system: the jet from the inner disc of the 
primary, and the molecular outflow from the inner bound- 
ary of the CB disc. The dust component of the azimuthally 
non-homogeneous disc wind could be a source of the vari- 
able extinction both in the case of HH 30, and other young 
stellar object^]. 
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2 Cotera et al. (2007) studying the images of the young stellar 
objects obtained with HST, found three other objects with the 
variable lateral asymmetries analogous to that detected in HH 30. 



